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Starting from the present version of the Riga dynamo experiment with its rotating magnetic eigenfield 
dominated by a single frequency we ask for those modifications of this set-up that would allow for a non- 
trivial magnetic field behaviour in the saturation regime. Assuming an increased ratio of azimuthal to axial 
flow velocity, we obtain energy oscillations with a frequency below the eigenfrequency of the magnetic field. 
These new oscillations are identified as magneto-inertial waves that result from a slight imbalance of Lorentz 
and inertial forces. Increasing the azimuthal velocity further, or increasing the total magnetic Reynolds 
number, we find transitions to a chaotic behaviour of the dynamo. 
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1 Introduction 

The last decade has seen great advantages in the exper- 
imental realization of hydromagnetic dynamos (Stefani, 
Gailitis and Gerbeth 2008), culminating perhaps in the 
recent observation of reversing and chaotic dynamos in 
the French VKS experiment (Ravelet 2008). 

The Riga dynamo experiment traces back to one of 
the simplest homogeneous dynamo concepts that had 
been proposed by Ponomarenko (1973). The Ponomer- 
anko dynamo consists of a conducting rigid rod that 
spirals through a medium of the same conductivity that 
extends infinitely in radial and axial direction. 

By a detailed analysis of this dynamo, Gailitis and 
Freibergs (1976) were able to identify the rather low 
critical magnetic Reynolds number of 17.7 for the con- 
vective instability. An essential step towards the later 
experimental realization was the addition of a straight 
back-flow concentric to the inner helical flow, which 
converts the convective instability into an absolute one 
(Gailitis and Freibergs 1980). 

After many years of optimization, design and and 
construction (cf. Gailitis et al. 2008 for a recent survey) 
a first experimental campaign took place in November 
1999, resulting in the observation of a slowly growing 
magnetic eigenfield (Gailitis et al. 2000). 

The four dynamo runs of July 2000 provided a first 
stock of growth rate, frequency, and spatial structure 
data of the magnetic eigenfield, this time both in the 
kinematic as well as in the saturated regime (Gailitis 



Corresponding author: e-mail: F.Stefani@fzd.de 



et al. 2001). During the following experimental cam- 
paigns, plenty of growth rate and frequency measure- 
ments were carried out. When scaled by the temper- 
ature dependent conductivity of sodium, both quan- 
tities turned out to be reproducible over the years. 
(Note, however, that after the replacement of the out- 
worn gliding ring seal by a modern magnetic coupler 
in 2007 a problem with a reduced azimuthal velocity 
occurred which is detrimental for dynamo action and 
which needs further inspection). 

With view on the future utilization of the Riga dy- 
namo experiment, here we deal with the general ques- 
tion whether the facility could be modified in such a 
way that it exhibits a more complicated saturation be- 
haviour. One idea that comes to mind is the princi- 
ple possibility to lower the critical Rm by means of 
a parametric resonance, or swing excitation, (Rohde, 
Riidiger and Elstner 1999) with a periodic flow forc- 
ing with approximately the double of eigenfrequency. 
We feel, however, that the comparably high frequency 
(around 4 Hz and larger) that would be needed for this 
purpose is far beyond the technical specifications of the 
utilized motors and their control. 

Another idea, namely to modify the present ratio 
of azimuthal velocity v v and axial velocity v z , is mo- 
tivated by the following considerations: The Riga dy- 
namo had been constructed in such a way that the 
flow helicity is maximum for a given kinetic energy 
(Gailitis et al. 2004) which implies the ratio of v v to 
v z to be close to 1. Now, since the saturation mech- 
anism relies strongly on a selective braking of v v one 
could ask for the consequences if one would start with 
a too strong v v . Its braking would then lead to a re- 
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Fig. 1 The Riga dynamo experiment and its eigen- 
field. (a) Details of the central dynamo module. 1 - 
Upper bending region; 2 - Propeller; 3 - Central helical 
flow region; 4 - Return-flow region; 5 - Outer sodium 
region; 6 - Guiding vanes for straightening the flow in 
the return flow. 7 - Lower bending region. At approxi- 
mately 1/3 (L) and 2/3 (U) of the dynamo height there 
are four ports for various magnetic field, pressure and 
velocity probes, (b) Simulated structure of the mag- 
netic eigenfield in the kinematic regime. 

duced ratio v v to v z which might be suspected to be 
even better suited for dynamo action than the original 
kinematic flow. In the extreme form this mechanism 
could lead to a sub-critical Hopf bifurcation, meaning 
that the velocity in the saturated state has a lower 
critical Rm than the kinematic (unperturbed) veloc- 
ity. This effect should be detectable in form of a hys- 
teretic behaviour of self-excitation, similar to the case 
studied recently by Reuter et al. (2008). But even if 
an " honest" sub-critical Hopf bifurcation could not be 
achieved, a general tendency of the system to oscillate 
between a weak-field and a strong-field state is still to 
be expected. Since such an oscillation would appear 
from the destabilization of the (former) steady equi- 
librium between Lorentz forces and inertial forces, it 
should be called a magneto-inertial wave, similar to the 
waves that have been observed in the DTS-experiment 
in Grenoble (Schmitt 2010). 

Further, if such an oscillation with a second fre- 
quency sets in, one could also ask for the appearance 
of more frequencies and possibly for a transition to 
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Fig. 2 Assumed kinematic velocity structure in the 
central part of the inner tube and in the back-flow tube 
at a rotation rate of 2100 rpm. Axial velocity v z (r, z) 
(l.h.s) and azimuthal velocity v v (r, z) (r.h.s.). Note the 
slight downward decay of v v (r, z). Figure adopted from 
(Gailitis et al. 2002) 



chaos according to the Ruelle-Takens-Newhouse sce- 
nario (Eckmann 1981). 

These foregoing ideas delineate the scope of the 
present paper. In the second section, we will describe 
briefly the Riga dynamo, and present the utilized scheme 
for its numerical simulation. Then we will discuss two 
sorts of parameter variations that both start from the 
well-known single-frequency saturation state and end 
up with a chaotic regime. The paper concludes with a 
few remarks concerning the technical feasibility of the 
proposed modification. 

2 The Riga dynamo and the numerical 
code for its simulation 

An illustration of the Riga dynamo, and of its mag- 
netic eigenfield is shown in Fig. 1. The facility con- 
sists basically of three concentric cylinders. In the in- 
nermost cylinder, the liquid sodium undergoes a down- 
ward spiral motion with optimized helicity. In the sec- 
ond cylinder, it flows back to the propeller region. The 
third cylinder serves for improving the electromagnetic 
boundary conditions which helps in achieving a mini- 
mum critical Rm. More details about the facility and 
the experimental campaigns can be found in (Gaili- 
tis et al. 2000, 2001, 2003, 2004, 2008). In the present 
configuration, self-excitation occurs approximately at 
a propeller rotation rate of P = 1840 rpm when scaled 
appropriately to a reference temperature of 157° C. The 
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Fig. 4 The three lowest dominant frequencies for 
increasing factor A. f\ denotes the frequency of the 
magnetic eigenficld, / 2 is the frequency of the first 
magneto-inertial oscillation, fs is the third arising fre- 
quency, fi corresponds approximately to the combina- 
tion /i + f 2 - h- 

shape of the magnetic eigenfield resembles a double- 
helix that rotates around the vertical axis with a fre- 
quency in the range of 1-2 Hz. 

The numerical simulations of this paper will be car- 
ried out by means of a code that couples a solver for 
the induction equation 

OH 1 

— =Vx(vxB) + AB . (1) 

at noa 

for the azimuthal m = 1 component of the magnetic 
field B, and a simplified solver for the velocity per- 
turbation. The solver for the induction equation is a 
two-dimensional finite difference scheme with a non- 
homogeneous grid in r and z, that uses an Adams- 
Bashforth method of second order for the time integra- 
tion (cf. Gailitis et al. 2004 for more details). The real 
geometry of the dynamo module has been slightly sim- 
plified, with all curved parts in the bending regions be- 
ing replaced by rectangular ones. The unperturbed ve- 
locity in the central channel was inferred from a number 
of measurements at a water test facility and some ex- 
trapolations. These measurements had revealed a slight 
downstream decay of v v . The velocity in the back- flow 
channel has been considered as purely axial and con- 
stant in r and z, with the constraint that the volume 
flux there is equal to the volume flux in the inner cylin- 
der. Figure 2 shows the assumed kinematic velocity 
components v z (r,z) and v v (r,z) both in the central 
channel and in the back-flow channel. 

The induction equation solver has been extensively 
used in the optimization of the Riga dynamo experi- 
ment, and its main results (in particular concerning the 
growthrate and the frequency of the eigenfield) agree 
to a very few percent with the experimental results. 

As for the saturation regime, our code tries to take 
into account the most important back-reaction effect 



within a simple one-dimensional model (see Gailitis et 
al. 2002, 2008). While v z can be assumed rather con- 
stant from top to bottom due to mass conservation, v v 
can be easily braked by the Lorentz forces without any 
significant pressure increase. In the inviscid approxi- 
mation, and considering only the m = mode of the 
Lorentz force, we end up with the differential equation 
for the Lorentz force induced perturbation 5v v {r,z) of 
the azimuthal velocity component: 
d d 

—5v v (r,z,t) = -v z (r,z)—5v ip (r,z,t) (2) 

+ — [(V x B) x B]Jr,z,t) . 
Hop 

In contrast to the procedure described in (Gailitis et al. 
2002, 2004, 2008) which was focused exclusively on a 
stationary saturated regime, here we take into account 
the time dependence of 5v v {r, z,t) explicitly. 

Note that Eq. 2 is solved both in the innermost 
channel where it describes the downward braking of 
v v , as well as in the back- flow channel where it de- 
scribes the upward acceleration of v v . Both effects to- 
gether lead to a reduction of the differential rotation 
and usually to a deterioration of the dynamo capabil- 
ity of the flow, quite in accordance with Lenz's rule. 
As a technical remark, note that Eq. 2 is time-stepped 
independently in the innermost channel and the back- 
flow channel, with time-independent profiles set on two 
disconnection lines. The first disconnection line is in the 
central channel directly behind the propeller where the 
Vtp(r) profile is set to the water-test profile (multiplied 
by A). The second disconnection line is at the bottom of 
the back-flow channel behind the straightening blades 
where v v (r) is always set to zero. 

The main result of this code, the downward accu- 
mulating braking of v v , has been shown to be in rea- 
sonable agreement with experimental data, and with 
the results of a much more elaborated T-RANS simu- 
lation (see Kenjeres et al. 2006, 2007), apart from the 
fact that the latter simulation results in an additional 
slight radial redistribution of v z that cannot be covered 
by our simple back-reaction model. 

In the following we will use the physical units of the 
Riga dynamo experiment, i.e. the real geometry, veloc- 
ity and conductivity with the only compromise that we 
will always use the conductivity a = 8.56 x 10 6 S/m 
at the reference temperature T — 157°C (which was 
used in all previous publications on the Riga dynamo) . 
If we would dare, however, to lower the temperature 
to T = 102°C (i.e. 5°C above melting temperature) 
we would get a — 10.28 x 10 6 S/m, and the presently 
achievable propeller rotation rate of P = 2500 rpm 
would than be equivalent (in terms of the magnetic 
Reynolds number) to 3000 rpm. 

With the reference conductivity a = 8.56 x 10 6 S/m 
we obtain a diffusion time tdis = = 0.168 sec if 

the innermost radius R\ = 0.125 is used, or 1.72 sec if 
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Fig. 3 Influence of the increasing ratio of v v and v z (factor A compared with the present flow structure) on 
the time evolution. From top to bottom: A=l, 1.3, 1.7, 2.0, 2.5, 3.0. From left to right: Magnetic field energy and 
B z (r = 0.086 m, z = —0.39 m) (first column), two-dimensional phase portrait of B z and v v at (r = 0.086 m, z = 
—0.39 m) (second column), three-dimensional phase portrait of B r , B v , and B z (r — 0.086 m, z = —0.39 m) 
(third column), power spectral density of B z (r — 0.086 m, z = —0.39 m) (fourth column). 



we take the outer radius i?3 — 0.4 m. The integration 
time for our runs was set to 6.25 sec, which in some 
cases is not completely sufficient for all relaxations to 
be finalized. 



3 Numerical results 

In this section we present the numerical results for two 
types of parameter variations. The first variation starts 
from a state that corresponds approximately to the 
maximum achievable Rm number of the present Riga 
dynamo set-up, from where we artificially scale up v v 
by a factor A, keeping v z fixed. In the second variation 
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Fig. 5 Same as Fig. 3, but for increasing propeller rotation rate P at fixed factor A = 1.5. From to to bottom: 
P = 1900, 2600, 3000, 3600, 4500, 6000 rpm. 



we scale-up the magnetic Reynolds number (in terms 
of the propeller rotation rate P) while keeping the ratio 
of v v and v z fixed at A = 1.5. 

The variation documented in Fig. 3 starts from a 
given propeller rotation rate P = 3000 rpm with the 
present ratio of v v to v z (corresponding to A = 1). 
Keeping the profile of v z fixed, we artificially increase 
then v v . The rows from top to bottom in Fig. 3 cor- 
respond to increasing factors A =1.0, 1.3, 1.7, 2.0, 2.5, 
3.0. For each row, the first column shows the time evo- 
lution of the magnetic field component B z at radius 



r = 0.086 m and axial position z — —0.39 m (with 
respect to the mid-height of the dynamo) and of the 
total magnetic field energy. The second column depicts 
the two-dimensional phase portrait of B z and v v at 
the same position (r = 0.086 m, z = —0.39 m). The 
third column shows a three-dimensional phase portrait 
in the three magnetic field components B r , B Vl and 
B z at the same position. Finally, the fourth column 
gives the power spectral density of B z (r = 0.086 m,z = 
-0.39 m) 
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oscillation (resulting from the rotating magnetic eigen- 
field) with a frequency /i = 2.13 Hz, which trans- 
lates into a (more or less) horizontal line in the B z - 
v v plot and a circle in the B r -B v -B z phase portrait. 
For A = 1.3, a second frequency /2 = 1.71 Hz comes 
into play, which is a signature of an arising magneto- 
inertial oscillation. Accordingly, the three-dimensional 
phase portrait shows a torus structure. Around A = 2.0 
a third frequency arises. For A = 2.5 and A = 3.0 we 
observe then a irregular, probably chaotic behaviour, 
which manifests itself in the phase portraits and the 
smeared out PSD. 

In Fig. 4 we have compiled the three lowest peaks 
/i, /2, and fs of the frequency spectrum. We have also 
indicated a fourth frequency, which seems to agree 
with the combination fi+fi—fs- It is important to note 
that these spectral properties have still a preliminary 
character because the simulation time comprises only 
a few oscillations. Anyway, there is some evidence that 
the appearance of the third frequency "rings in" the 
on-set of irregular behaviour by destabilizing the two- 
torus made up by the first two frequencies /i and fa- 

A second parameter variation, for which the value 
A is now fixed to 1.5, while the propeller rotation rate 
P is increased, is illustrated in Figs. 5 and 6. The rows 
from top to bottom in Fig. 4 correspond now to the ax- 
ial propeller rotation rates P =1900, 2600, 3000, 3600, 
4500, 5000, 6000 rpm. 

The main features are quite similar as in the former 
sequence, apart from the fact that the ratio /1//2 ~ 2 
stays rather constant within the wide range of P = 
1900. ..4000 rpm, a fact that points to a phenomenon 
called "mode locking". This is best visible in the ap- 
pearance of a rather sharp ratio /1//2 = 2.0 for P = 
3600 rpm (fourth row) in Fig. 5). 

Figures 7 and 8 illustrate the character of the magneto- tions of the sodium flow in which the present ratio of 
inertial wave for A = 1.5 and two selected values P = v v and v z is scaled-up by a factor A. Beginning ap- 
3600 and 6000 rpm, respectively. In each case, the up- proximately at A = 1.3 we observe the appearance of 
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(b)0 0.5 t 1 1.5 

Fig. 7 t — z-dcpcndence of v v (r = 0.086 m, z, t) (a) 
and B r (r = 0.086 m, z, t) (b) for A = 1.5 and P = 3600 
rpm. The unit of v v is m/s, the unit of B r is Tesla. 

per panel shows v v {r = 0.086 m, z,t), and the lower 
panel shows B r (r = 0.086 m, z, t). While the case P = 
3600 rpm shows a frequency-locked wave, the case P = 
6000 is already much more irregular. 

4 Conclusions 

Presently, the maximum Rm in the Riga dynamo ex- 
periment can be achieved with a propeller rotation rate 
of 2500 rpm and a temperature of around 102° C, which 
is roughly equivalent to a propeller rotation rate of 3000 
rpm at a reference temperature of 157°C. Starting from 
this realistic set-up that leads to a single- frequency sat- 
urated state, we have studied hypothetical modifica- 



© 2011 WILEY- VCH Verlag GmbH & Co. KGaA, Wcinhoim 



www.an-journal.org 



Astron. Nachr. / AN (2011) 



7 



VP : -36 -34 -32 -30 -28 -26 -24 -22 -20 -1 8 -1 6 -1 4 -1 2 




BR Real: -o.os -0.06 -0.04 -0.02 0.02 0.04 0.06 o.os 




(b)0 0.5 t 1 1.5 

Fig. 8 Same as Fig. 7, but for A = 1.5 and P = 6000 
rpm. 

a magnetic field energy oscillation with a frequency / 2 
that is in the order of one half of the eigenfrequency /1 
of the saturated dynamo. 

We hypothesize that this new oscillation with fre- 
quency /2 is basically a magneto-inertial wave that 
traces back to the fact that the steady equilibrium be- 
tween Lorentz and inertial forces (that is responsible 
for the saturation) becomes unstable if the ratio of v v 
and v z velocity is larger than the optimum value. This 
allows the system to oscillate between a weak field and 
a strong field limit. 

For increasing factors A we observe further the ap- 
pearance of a third frequency fz and then the transition 
to a rather irregular, probably a chaotic state. Similar 
transitions are observed when increasing the propeller 
rotation rate P while holding the factor A fixed at a 
value of 1.5. An interesting feature of this latter se- 
quence is the appearance of a sort of frequency locking 



(characterized /1//2 = 2 ) for a rather wide range of 
2600 < P < 3600. 

The observed behaviour resembles the Ruelle-Takcns- 
Newhouse scenario of the transition to chaos (see Eck- 
mann 1981). We have to point out, though, the pre- 
liminary character of our study, both with respect to 
limited integration time as well as to the fact that a 
detailed study of the space and time resolution is still 
missing. 

What remains to be checked further is the possible 
existence of a sub-critical Hopf-bifurcation. 

The delineated mechanism to obtain energy oscil- 
lations in the saturated state, and even the sketched 
route to chaos, seems to be in the range of technical 
feasibility if we assume approximately a doubling of the 
motor power and an increase of the ratio of azimuthal 
to axial velocity by means of an appropriate propeller 
and guiding blade design. 
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